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FINAL REPORT 


ONR “OCEAN SCIENCE EDUCATOR IN NUMERICAL MODELING" 

Four projects, all leading to publication, were supported by the ONR Educator 
Award. 

1) The first is an investigation of Gulf Stream dynamics using concepts of chaotic 
transport. A Gulf Stream-like meandering current is kinematically modeled using a 
jet-like streamfunction. The streamfunction has been used by other authors to 
investigate geophysical fluid motion. Under suitable parameter values a "front" 
exists along the velocity maximum across which there is no transport, but with 
mixing occurring on either side. An analogous barrier exists in the Gulf Stream, 
but appears to fade with depth, as shown by others. The conditions for mixing 
across the model front are studied using the Chiiikov oveiiap criterion which 
indicates the front breaks down for meander amplitude above a critical threshold 
that is inversely dependent on the ratio of meander phase speed and current speed, 
c/U. 

It appears that the increase in cross-frontal mixing in the deeper levels of the 
Gulf Stream is the x'esult of current meandering and the deci'ease of current velocity 
with depth. The mechanism is the interaction of difiez-ent meander modes traveling 
along the Gulf Stream. 

2) The second topic was a discussion of the sensitivity of objective techniques to 
parametric variations. Solutions of some objective analysis techniques are known to 
depend upon the subjective values of internal parameters. The change in the 
solution per change in the parameter is the sensitivity. Parameters with low 
sensitivity can be varied with laz'ge increments during preliminary searches for 








near-optimal parameter values. Only terms with high sensitivity must be 
thoroughly investigated once the parameters are determined to be close to optimal. 
Both absolute and relative sensitivities are discussed and a sensitivity-based 
definition of the solution uncertainty is proposed. The sensitivity of direct 
minimization analysis to parametric variation is evaluated using a set of "response 
function" that characterize different aspects of the solution. 

Two examples are used to illustrate the usefiilness of the technique. Both 
involve measurements of air-sea quantities (e.g., wind stress and latent heat flux) 
from a variety of data sources using a direct minimization technique. The examples 
demonstrate that sensitivity analysis is capable of quantifying regional sensitivities 
as well as indicating the magnitude and relationship between the various 
parameters. 

3) The third topic is an early application of wavelet methods to equatorial waves. 
Wavelet analysis is a relatively new technique that is an important addition to 
standard signal analysis methods. Unlike Fourier analysis which yields an average 
amplitude and phase for each harmonic in a data set, the wavelet transform 
produces an "instantaneous" estimate or local value for the amplitude and phase of 
each harmonic. This allows detailed study of nonstationary spatial or time- 
dependent signal characteristics. 

We demonsti-ate the usefulness of the transform by studying the dispersion of 
Yanai waves in a reduced gravity equatorial model. The gr oup velocity is measured 
directly over a finite range of wavenumbers by examining the time evolution of the 
transform. The results agree well with linear theory at higher wavenumber but the 
measured group velocity is reduced at lower wavenumbers, possibly due to 
interaction with the basin boundaries. 










4) The fourth study is an application of wavelet analysis. Two-year sea-surface 
temperature time series of satellite data at two sites in the equatorial Indian Ocean 
are examined for oscillations with periods of 2-70 days. The wavelet transform of 
the signals reveals a changing wavelet spectrum between August 1987 and 
November 1987 in the 10-30 day range, while the same time period in 1988 shows a 
relatively fixed spectrum. At 3 degrees latitude and 53E longitude the 1987 wavelet 
coefficients with scales 10-30 days have about twice the amplitude they have in 
1988. At 3 latitude and 56E longitude the 1987 waves have roughly half the 
amplitude as the 1988 waves. Activity with wavelet spectral peaks at periods near 
12 says often precede these waves. 
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